Spontaneous Breakdown of Translational Symmetry in Quantum Hall Systems: 

Crystalline Order in High Landau Levels 
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We report on results of systematic numerical studies of two-dimensional electron gas systems subject 
to a perpendicular magnetic field, with a high Landau level partially filled by electrons. Our results 
are strongly suggestive of a breakdown of translational symmetry and the presence of crystalline 
order in the ground state. This is in sharp contrast with the physics of the lowest and first excited 
Landau levels, and in good qualitative agreement with earlier Hartree-Fock studies. Experimental 
implications of our results are discussed. 
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Recently there has been considerable interest in the be- 
havior of a two-dimensional (2D) electron gas subject to 
a perpendicular magnetic field, when a high Landau level 
(LL) (with LL index > 2) is partially filled by elec- 
trons. This is largely inspired by the recent experimental 
discovery ||-||] that the transport properties of the sys- 
tem are highly anisotropic and non-linear for LL filling 
fraction v = 9/2, 11/2, 13/2, •• ■. Previously, Hartree- 
Fock (HF) and variational studies Q suggested 
that, unlike the = and N = 1 LL's (in which either 
incompressible fractional quantum Hall (FQH) or com- 
pressible Fermi-liquid like states are realized), in iV > 2 
LL's the electrons form charge density waves (CDW). In 
particular, at half-integral filling CDW's break transla- 
tional symmetry only in one direction and form stripes. 
Anisotropic transport would indeed result from such a 
striped (or related) structure |P-pT|]. 

We neglect LL mixing, and consider the case where 
the LL with index N has partial filling i>, while LL's 
with lower index are completely filled (y — 2N + D). By 
particle-hole symmetry of the partially-filled LL, this is 
equivalent to = 2A^ + 2 — D. We also assume that 
the partially-filled LL is maximally spin-polarized at the 
i> we consider. Previously |lj], we studied such N > 
2 LL's with = 1/2 by numerically diagonalizing the 
Hamiltonians of finite-size systems; those results strongly 
supported the existence of stripe order. 

An outstanding issue is the nature of the ground state 
at high LL's for fillings sufficiently far from the half-filled 
level. Koulakov, Fogler and Shklovskii |^,|| (see also 
Moessner and Chalker Q) predicted a novel crystalline 
phase called the "bubble" phase with more than one elec- 
tron per unit cell outside of the range p = 0.4 — 0.6. The 
bubble crystal has lower energy than the Laughlin state 
for = 4+1/3 01 . Experimentally, a re-entrant quantum 
Hall state is found near = 4 -|- 1/4 which is quantized 
as a = 4 LL plateau Evidently the electrons in 

the top- most LL are frozen out of the transport. Pinning 
of a crystalline structure provides a natural explanation 



of the re-entrant phase and would further explain the ob- 
served threshold in conduction However this is not 
entirely conclusive and other mechanisms for the conduc- 
tion threshold are also possible O]. 
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FIG. 1. Energy levels versus aspect ratio for quarter-filled 
TV = 2 Landau level with eight electrons and rectangular ge- 
ometry. The inset is a blow-up of the low-energy spectra for 
aspect ratio between 0.6 and 1.0. The points at 1.1 (open 
circle plus x) correspond to hexagonal unit cell. 

In this paper, we report on new numerical results 
on systems away from half-filling using the unscreened 
Coulomb interactions. Remarkably, our results suggest 
that CDW's are formed at all filling factors we have 
studied, including those that would support prominent 
FQH states or composite fermion Fermi-liquid states in 
the lowest or first excited Landau levels. These CDW's, 
however, have 2D structures and are no longer stripes 
when the filling factors are sufficiently far away from 1/2. 
They are not Wigner crystals either, unless D is small 
(below 0.2). In the intermediate filling factor range, we 
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find each unit cell of the CDW contains more than one 
electron. Our results are in good agreement with the 
predicted bubble phase and are the first exact finite-size 
calculations which exhibit a crystalline state in a system 
with continuous translational symmetry. 
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FIG. 2. The allowed (circles) and the ground state mani- 
fold momenta (x's). The data is for a rectangular geometry 
with asp = 0.77, 8 electrons in the N = 2 LL and v = 1/4. 
The solid line is the boundary of the BZ. The superlattice 
reciprocal basis vectors are shown by solid arrows. The inset 
gives the corresponding results for an hexagonal unit cell. 

We restrict the states of the electrons to a given LL, 
and work with periodic boundary conditions (PBC, torus 
geometry) as in our previous paper We also set 

the magnetic length to unity. To detect intrinsically pre- 
ferred configurations we consider a rectangular PBC unit 
cell and vary its aspect ratio. The PBC plays a crucial 
role in removing continuous rotational symmetry, and se- 
lecting a discrete set of possible crystal orientations. 

In Fig.0 we plot the energy levels of systems with 
iVe = 8 electrons in the N = 2 LL at filling factor i> — 1/A 
as a function of the aspect ratio. We also show the levels 
of a system with a hexagonal PBC unit cell at the right 
side of Fig.^ A generic feature of the spectra is the exis- 
tence of a large number of low- lying states whose energies 
are almost degenerate, which we call the ground state 
manifold. The momenta of these quasi-degenerate states 
for rectangular geometry with aspect ratio asp = 0.77 
and hexagonal geometry are shown in Fig.||; they form a 
2D superlattice structure, which for the rectangular ge- 
ometry have the super cell vectors bi — 2aex — bSy, and 
b2 — 2bey. Where, a ~ 2tt/Li, and b = 2tt/L2. Li, and 
L2 are the dimensions of the unit cell (Li x L2 ~ 2nN^, 
is the total flux quanta in the system) . The area per 
wavevector in the Brillouin zone (BZ) is ab = (27r)^/A, 
where A is the (real space) area of the system. 



There are similarities as well as important differences 
between these spectra and those of half-filled high LL's 
]l2t with stripe order. As in the stripe case the large 
quasi-degeneracy of the ground state manifold is an in- 
dication of broken translational symmetry ]l5[ | . The dif- 
ference here is that (i) the degeneracy is much larger and 
(ii) the momenta of the low-lying states form 2D instead 
of ID arrays. These new features indicate that the trans- 
lational symmetry is broken in both directions and the 
ground state is a 2D CDW. In the stripe state, on the 
other hand, the translational symmetry is only broken in 
the direction perpendicular to the stripes. Therefore the 
degeneracy is smaller and the momenta of the low-lying 
states form a ID array. 
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FIG. 3. A 2D plot of the peaks of the projected (or guiding 
center) charge susceptibility x(ci) Eit reciprocal lattice vectors, 
for a system with Ne — 8 ai v — 1/4, rectangular geometry 
with a = 0.77. The size of the circles give an indication of the 
height of the peak at that point. Only responses above 100 
have been plotted as solid circles. The zone boundaries are not 
within the range of the figure. The largest circle corresponds 
to 16491. The inset gives the results for a hexagonal unit cell. 
The x's are the allowed wavevectors. 

The momenta of the states in the ground state mani- 
fold are the reciprocal lattice vectors of the bubble crys- 
tal. Transforming to the direct lattice vectors, we obtain 
ai — n/aex and a.2 ~ Tr/2aex + n/bSy. For the optimum 
system, with asp— 0.77, we obtain oi — 8.08, 02 — 7.42, 
and (f) — 57°. This is very close to a triangular lattice. In 
the case of hexagonal PBC unit cell, both the reciprocal 
superlattice and its direct lattice are triangular. 

The number Nd of distinct quasi-degenerate ground 
states allows the number Nh of bubbles in the system, and 
hence the number M = N^/Nb of electrons per bubble, 
to be immediately obtained through the relation NbNo 
— , where N is the highest common divisor of and 
iV$. In our case, N = = 8, and Nd = 16, which gives 
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Nf, = 4: and M = 2 . The Wigner crystal would corre- 
spond to Nb = Ne and M = 1. In general |17|, there are 
N'^ distinct values of the total momentum quantum num- 
ber, which define a BZ of area (27riV)^/A. If translational 
symmetry is broken, the area of the BZ of the superlat- 
tice is then {2TrN)'^ /ANd, which must be (27r)V(^/^(,), 
where A/Ni, is the area per bubble; hence iVf,A^£) = iV^. 



also seen in the "guiding center (GC) static structure fac- 
tor" S'o(q) 11, which we present as a 3D plot in FigJJ. 
Here we see sharp peaks with an approximate six-fold 
symmetry at the primary reciprocal lattice vectors, indi- 
cating the presence of strong density correlation at these 
wave vectors in the ground state. 




FIG. 4. A 3D plot of the guiding center structure factor 
So(q) (same system as Fig.^. The signature of the hexagonal 
lattice is seen in the near six-fold symmetry of the peaks. 

We next turn to the density response functions. In 
Fig.^ we show the projected ground state charge sus- 
ceptibility of one of the optimum rectangular 
and hexagonal systems described above. The calculation 
takes into account the contributions from the two lowest 
energy states in each symmetry subspace; this is an excel- 
lent approximation in view of the fact that the response 
function is dominated by low energy states because of the 
energy denominator. We note that x(q) exhibits a strong 
response at the reciprocal lattice vectors (Bragg condi- 
tion); the background at other wave vectors (shown by 
X 's in Fig.^) are negligible compared to these responses. 
The origin of the strong response lies in the approximate 
degeneracy among the states forming the ground state 
manifold. The system responds very strongly to a po- 
tential modulation with a wave vector that connects the 
ground state to one of the low lying states (which must 
be a reciprocal lattice vector) because of the small en- 
ergy denominator. This is also another reason why there 
must be one low-lying state for each reciprocal lattice 
vector. A second notable feature is the almost hexagonal 
symmetry of the response, despite the fact that the PBC 
geometry used in this case was rectangular. This indi- 
cates that the bubbles tend to form a triangular lattice, 
in agreement with the predictions of HF theory. 

The tendency toward forming a triangular lattice is 




FIG. 5. Real-space "projected density" (guiding center) 
correlation function, derived from Fig.^. A 2D contour plot 
is also included below the 3D plot. 




FIG. 6. A'^ = 2 Landau level full electron density correla- 
tions (relative to a guiding center), derived from Fig.^. 

In Fig.||and FigJ| we plot ground state "projected den- 
sity" correlation functions in real space. These describe 
correlations relative to the GC {not the coordinate) of 
a particle. The first is the Fourier transform (FT) of 
S'o(q) exp(— g^/2), which is the electron density of an 
equivalent lowest-LL system, and gives information on 
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the spatial distribution of GC's. The second is the FT of 
S'o(q)[Ljv(gV2)]^exp(-gV2) with N = 2 {Lm is a La- 
guerre polynomial): this (plus the uniform density of the 
filled LL's) represents the actual electron density. 

In Fig.| the presence of four bubbles and the relative 
orientation of the bubbles can be clearly seen and there 
is strong crystalline order of the GC distribution. The 
central peak contains two electrons, one of which is the 
particle with the GC at the origin. For > 0, as in 
Fig.^, only weak order is displayed by the actual electron 
density, because of the averaging effect of the cyclotron 
motion around the GC's. It is the guiding centers of the 
electrons that form bubbles as anticipated in reference 7 
(Fig.|l|). The electrons themselves manage to stay apart 
to lower the Coulomb repulsion, in spite of the clustering 
of their GC's. 




always pinned by a disorder potential, and in a nonlinear 
transport measurement, there should be a threshold de- 
pinning field at which there is a sharp feature in the I — V 
curve. A weaker crystal would result in more diffuse con- 
duction threshold as various portions of the crystal get 
depinned at different current values, while a stronger one, 
on the other hand, will have sharp conduction threshold. 
This is consistent with the observation of Cooper et al |^ 
that there is a sharp threshold region at about — 1/4, 
but more diffuse thresholds at both higher and lower v. 

In contrast to i/ = 1/4 and i' = 1/3, the spectra for 
v — 2/b was found to be very similar to = 1/2 |l^]. The 
momenta of the low-lying states belong to a ID array, 
indicating formation of a ID CDW or stripe phase; the 
weight of the HF state 11110000001111000000 {N^ = 8), 
in a rectangular geometry with an aspect ratio of 0.80, is 
about 65%. We conclude that the transition from stripe 
to bubble phases occurs between i> = 1/3 and i' — 2/5, 
in qualitative agreement with HF predictions [§j6|. We 
have also studied higher LL's. The results are similar 
and will be reported elsewhere. 
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and DMR-0086191 (E.H.R.), DMR-9809483 (F.D.M.H.), 
DMR-9971541, and the Sloan Foundation (K.Y.). E.H.R. 
acknowledges the hospitality of ITP Santa Barbara, sup- 
ported by NSF-PHY94-07194, where part of the work 
was performed. 



FIG. 7. Spectra of systems with eight electrons at i/ = 1/3 
in the N — 2 LL, with rectangular geometry and various 
aspect ratios. The inset plots the momenta of the low-lying 
states for asp — 0.75, in the same way as Fig. ^ 

We have also explored other filling factors in the N = 2 
LL such as = 2/5 and z> = 1/3, where the system would 
condense into prominent FQH states if it was in the low- 
est LL; here, however, our studies suggest formation of 
CDW's instead. For i/ = 1/3 we obtain similar behav- 
ior to = 1/4: the energy spectra as a function of the 
aspect ratio shown in Fig.^ is very similar to Fig.|^ and 
indicates formation of a 2D CDW. Using the degeneracy 
of the ground state manifold, we find the number of elec- 
trons per bubble is also two. The energies of the states 
in the ground state manifold, however, are not as close 
as the i> — 1/4. This results in weaker peaks in x(q) 
the reciprocal lattice vectors. 

We interpret this to be an indication that, in this LL, 
i> = 1/4 is more favorable than i> = 1/3 for formation of 
a two-electron bubble phase. In real systems a crystal is 



[1] M. P. Lilly, K. B. Cooper, J. P. Eisenstein, L. N. Pfeiffer, 
and K. W. West, Phys. Rev. Lett. 82, 394 (1999). 

[2] R. R. Du, D. C. Tsui, H. L. Stormer, L. N. Pfeiffer, K. 
W. Baldwin, and K. W. West, Solid State Commun. 109, 
389 (1999). 

[3] K. B. Cooper, M. P. Lilly, J. P. Eisenstein, L. N. Pfeiffer, 
and K. W. West, Phys. Rev. B 60, R11285 (1999). 

[4] A. A. Koulakov, M. M. Fogler, and B. I. Shklovskii, Phys. 
Rev. Lett. 76, 499 (1996). 

[5] M. M. Fogler, A. A. Koulakov, and B. I. Shklovskii, Phys. 
Rev. B 54, 1853 (1996). 

[6] R. Moessner and J. T. Chalker, Phys. Rev. B 54, 5006 

(1996) . 

[7] M. M. Fogler and A. A. Koulakov, Phys. Rev. B 55, 9326 

(1997) . 

[8] E. Fradkin and S. A. Kivelson, Phys. Rev. B 59, 8065 
(1999); E. Fradkin, S. A. Kivelson, E. Manousakis and 
K. Nho, Phys. Rev. Lett. 84, 1982 (2000). 

[9] H. A. Fertig, Phys. Rev. Lett. 82, 3693 (1999). 
[10] A. H. MacDonald and M. P. A. Fisher, Phys. Rev. B 61, 
5724 (2000). 

[11] T. Stanescu, I. Martin, and P. Phillips, Phys Rev. Lett. 



84, 1288 (2000). 
[12] E. H. Rczayi, F. D. M. Haldane, and K. Yang, Phys. Rev. 

Lett. 83, 1219 (1999). 
[13] H. Fukuyama, P. M. Platzman, and P. W. Anderson, 

Phys. Rev. B 19, 5211 (1979). 
[14] F. D. M. Haldane, in The Quantum Hall Effect, 2nd Ed., 

edited by R. E. Prange and S. M. Girvin (Springer, New 

York, 1990). 

[15] True broken translational symmetry with exact ground 
state degeneracy only occurs in the infinite-size hmit; at 
finite system size, this becomes a quasi-dogoncracy. 

[16] This result is in overall agreement with the estimate of 
reference 5: M = 3uN. 

[17] F. D. M. Haldane, Phys. Rev. Lett. 55, 1095 (1985). 



5 



